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Abstract 

Data on neutrino oscillation often involve reconstructed neutrino energies while the analysis 
implies the real neutrino energy. The corrections corresponding to the transformation from real to 
reconstructed energy are discussed in the case of Cherenkov detectors where multinucleon events 
appear as quasielastic ones. These corrections show up as a tendency for the events to escape 
the region of high flux, with a clear preference for the low energy side. This is an effect of the 
multinucleon component of the quasielastic cross section. We have applied our corrections to the 
T2K and MiniBooNE data for electron appearance or disappearance data. We show that the 
inclusion of this correction in the analysis is expected to lead to an increase of the best fit oscillation 
mass parameters, particularly pronounced for the MiniBooNE neutrino data. This inclusion in 
the analysis of the MiniBooNE neutrino data should improve the compatibility with the existing 
constraints. 

PACS numbers: 13.15.+g, 25.30.Pt, 14.60.Pq 
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I. INTRODUCTION 



The interpretation of neutrino oscillation experiments requires the determination of the 
neutrino energy which enters the expression of the oscillation length. This determination is 
often done through "quasielastic events" and the reconstructed energy hypothesis in which 
the neutrino interaction in the nuclear target with lepton production is supposed to take 
place on a nucleon at rest. In this case the only measurement of the lepton variables, its 
energy and emission angle, allows the reconstruction of the neutrino energy. In Cherenkov 
detectors "quasielastic events" are defined as those in which the emission product only in- 
clude one lepton, the ejected nucleons being unobservable for these types of detectors. In 
this case the definition of quasielastic events incorporates multinucleon ones which are indis- 
tinguishable from the genuine quasielastic ones. We have suggested 1] that the MiniBooNE 
axial mass anomaly [2] in the quasielastic cross section could be accounted for by the in- 
clusion of the multinucleon channel and we have been able to reproduce the total 
"quasielastic" cross section as well as the measured double differential cross section with- 

n 

out any modification of the axial mass. In a more recent work [4| we have addressed the 
question of the energy reconstruction, taking into account the fact that not all events being 
real quasielastic ones, the usual reconstruction method becomes questionable. We have ex- 
plored, for charged current neutrino reactions, the corrections to this method in the form of 
a probability distribution, F(E U , E u ), to have a real energy E v starting from a reconstructed 
value E v . We have shown that this distribution can be expressed in terms of the double 
differential neutrino nucleus cross section with respect to the energy u> transferred to the 
nucleus and the lepton emission angle 9, ^fj^ e . The distribution also involves the neutrino 
flux distribution $>(E V ). Similar approaches taking as well into account the multinucleon 
contribution have followed js,^]. 

The double differential cross section dE d2 ^ cose for muon production has been measured 
by MiniBooNE 2j|. Although the ponderations by the neutrino flux do not allow a direct 
insertion of this experimental data to derive the energy distributions, the fact that our theo- 
retical model is able to account in a satisfactory way for the MiniBooNE data of this double 
differential cross section [3] allows a certain degree of confidence in the distributions that 
we have obtained. The most spectacular influence of the multinucleon component in the 
distribution is the existence of a tail in the neutrino energy region above the reconstructed 
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energy. With the evaluation of these probability distributions we could transform the ex- 
perimental distributions such as those given by MiniBooNE HQ or T2K Q, 0, in terms 
of reconstructed neutrino energies, to smeared ones expressed in terms of the true energies 
which may then be confronted to oscillation models. We have applied jj] the smearing pro- 
cedure to the oscillation electron events issued from the oscillation of a muon neutrino beam 
in the case of the T2K experiment where electrons from the interaction of electron neutrinos 
attributed to an oscillation process of a muon neutrino beam have been observed and their 
distribution in terms of the reconstructed energy given [h]]. Here the disappearance effect 
is well established and the oscillation parameters known to a good accuracy. Hence the 
electron neutrino energy distribution which enters the smearing function is predictable from 
the muon neutrino one. We have applied our smearing procedure to extract the distribution 
in terms of the real neutrino energy. After this procedure the shape of the distribution could 
be successfully compared with the expected theoretical predicted shape, product of our total 
quasielastic cross section with the electron neutrino distribution generated by the oscillation 
of muon neutrinos. This consistency confirms the validity of the assumptions made: the 
observed electron events do come from an oscillation phenomenon and the oscillation mass 
parameter has the right magnitude. The second case is the MiniBooNE one where excess 
electron events attributed to the oscillation of the muon beam have also been observed Isj]- 
Here one deals with a short baseline experiments, L = 541 m, and the grounds for describing 
these oscillations are not as established. Short baseline oscillation phenomena, if they exist, 
involve sterile neutrinos with large values of the Am 2 , in the eV 2 range. But this value, 
hence the electron neutrinos energy distribution, is uncertain. We remind that our smearing 
procedure depends on this distribution. Assuming as an example that it followed the muon 
neutrino one we have applied the smearing procedure to the experimental MiniBooNE distri- 
bution. The data present a low energy anomaly with an excess of low reconstructed energy 
electron events. After our smearing procedure the excess of low energy electron neutrino 
events which raises a problem in the oscillation models, is pushed towards higher energy, 
making it possibly compatible with oscillation models. Of course this conclusion was not 
reliable as the assumption made on the electron neutrino spectrum was not based on oscilla- 
tion models and therefore totally arbitrary. Indeed in this case we observe that, contrary to 
the T2K case, there is no consistency at all between the smeared experimental distribution 
and the theoretical one, product of the total charged current "quasielastic" cross section for 
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electron neutrino by their energy distribution, a Ue (E Ue )^>(E Ue ), the theoretical one extending 
much further in energy. This is not a surprise since the flux model was totally arbitrary. 
However this study indicated the possibility of improving the description of the anomaly 
through the reconstructed energy corrections. In the present work we discuss in a more 
realistic way the influence of the smearing procedure on the two experiments both for the 
muon events from the beams and for the electron ones coming from the oscillation. The 
aim is to show the possible influence of nuclear physics effects on the determination of the 
oscillation parameter. For MiniBooNE we discuss the cases both neutrinos and antineutri- 
nos. In the first one we have the information that our description of the neutrino nucleus 
interaction is able to reproduce the double differential cross section with respect to the muon 
energy and emission angle. In the case of antineutrinos no such cross check can be made. 
The total "quasielastic" cross section is not published nor is the double differential cross 
section. We have to rely on a purely theoretical approach, which was described in our Refs. 
. Its main features are summarized in the following section. 

II. THEORY 

In the following the word neutrino will mean neutrino or antineutrino, the lepton can be 
either a fi~' + or an electron (positron), depending on the reaction in question. For Mini- 
BooNE the oscillation events concern electron neutrinos and the interesting events represent 
the production of electron arising from the interaction of electron neutrinos in the Cherenkov 
detector. The neutrino cross sections on nuclear targets are expressed in terms of nuclear 
responses, which represent the inelastic cross sections for a given type of coupling of the 
probe. For the charged current interaction case the probes have a purely isovector char- 
acter and they can involve or not the nucleonic spins. Our description of Refs. ll I3I. \l2\ 
introduces the multinucleon component only in the spin isospin response. For the com- 
parison between neutrinos and antineutrinos the important point is the role played by this 
isovector response relative to that of the spin isospin ones. Due to the negative sign of the 
axial vector interference term for antineutrinos the spin isospin contribution is diminished 
in the last case with respect to the isovector component. Consequently in our description 
the relative role of the multinucleon part is smaller for antineutrinos. This particular point 
of our treatment of the multinucleon channel has been challenged [l3, [yj]. Experimental 
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data will select between different approaches but, whatever the outcome, the method that 
we introduce for extracting a real distribution energy distribution from the reconstructed 
one and our general conclusions on the reconstruction effect remain valid. The difference 
of nuclear effects between neutrinos and antineutrinos is a very important point as it can 
mimic or mask the display of CP violating effect which introduces an asymmetry between 
v and v. 

In our last work j^J we started from the electron events experimental distributions in 
terms of the reconstructed energy. We transformed these into distributions in terms of the 
real neutrino energy. We remind that this transformation itself depends on the electron 
neutrino energy distribution hence on the oscillation parameter. Then we evaluated the 
theoretical expectations of the distribution with the same oscillation parameter. We could 
then confront theory and experiment. Although perfectly valid this method has the following 
drawback. When data evolve, the whole smearing procedure has to be redone. Moreover as 
each reconstructed energy bin influences the whole real energy distribution, negative values 
in certain bins influence the whole distribution. Here we will proceed differently. Indeed, as 
we observed in [4|, the procedure is completely reversible and can be used in both directions. 
Here we calculate the theoretical prediction for electron events energy distribution for a 
given value of the oscillation parameter. We then transform this distribution into one in 
terms of the reconstructed energy value, which can be directly compared to the experimental 
distribution. In principle we are then in a situation to investigate which oscillation parameter 
best fits the data. We also apply our smearing procedure to disappearance effects for the 
muon neutrinos in the T2K beam. 



A. Formalism 



The number of charged current events in a target for neutrinos of energy between E v and 
E v + dE u , for an energy transferred to the nuclear system, u, and a lepton emission angle 
9, is related to the double differential cross section by : 

d 2 a 

g(E u ,u,cos8) dE v du dcos6 = - — §(E U ) dE v du dcosO. (1) 

duj dcosO 

The quantity g is the triple density, in terms of the three variables, E u , u and cos#. For 
our problem it is convenient to switch to another set of variables, E v ,Ei{th.e energy of the 
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lepton produced) and the reconstructed neutrino energy E v . The relations between the two 
set of variables are firstly : u = E v — E\. In addition cos 9 is related to the new variables Ei 
and E v by: 



E v Pi cofiO + M(E V - Ei) - E V E X + 



mi 



0. 



(2) 



where Pi is the lepton momentum, mi the charged lepton mass and M the nucleon mass. 
The modulus of the Jacobian for these variables transformation is (MEi — mf /2)(E u Pi)~ 1 , 
and the new density G(E V , Ei, E v ) is 

d 2 a 



G(E V , E h E v ) dE v dE t dE v = dE u dE l dE u MEl m '/ 2 

eIp 



du dcos8 



u)=E v —Ei, cos0=cos 8{Ei,E v ) 

(3) 



where cos8(Ei, E u ) is the solution of Eq. (j2j). After integration over the lepton energy 
this density can be used in both directions. Either to extract a distribution in terms of 
the real neutrino energy from a distribution in reconstructed energies, as was done in our 
previous work J4] where we had used normalized probabilities. Or, in the opposite direction, 
we start from a theoretical distribution expressed with real energies then we perform the 
smearing procedure to deduce the corresponding distribution of the events in terms of the 
reconstructed energy. For this we integrate over the lepton energy and over the real neutrino 
energy distribution, which provides the distribution, D rec (E u ), in terms of the reconstructed 
energy which can be compared to the data 

d 2 a 



P rec^Pv 1 



dE v $(E v ) 



fpmin 



iE ME,- m y2 
EP, 



du dcos9 



(4) 



ui=E v —Ei, cos0=cos 9(Ei,E v ) 



where the quantities EJ nm and E™ ax are the minimum and maximum values of the charged 
lepton energy for a given value of E v . They are obtained by taking cos# = 1( — 1) in the 
Eq. 02]), with the additional restriction, mi < E x < E u . The second integral on the r.h.s. 
of Eq. (BJ, which represents the spreading function, depends on E v and E v \ we denote it 
as d(E u , E v \ We give in Fig. [T]some examples of its E v dependence for several E v values. 
The np-nh low energy tail is the counterpart, in these variables, of the high energy one that 
we stressed in our previous work {4]. 
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E v (GeV) 



FIG. 1: (Color online) The spreading function d{E v ,E v ) of Eq. per neutron of 12 C in the 
case of electrons evaluated for three E v values. The genuine quasielastic (dashed lines) and the 
multinucleon (dotted lines) contributions are also shown separately. 



III. APPLICATIONS 



T2K 



Here the situation is relatively simple as one deals with a long baseline experiment 
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with oscillation mass parameters already known to a good accuracy. We have pointed out 
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4[ the interest of the study for T2K of the muon events spectrum both in the close detector 
and in the far detector since the two corresponding muonic neutrino beams have different 
energy distributions. The study of the reconstruction influence on the electron events in 
the far SuperKamiokande detector was performed in our Ref. [4j, it is discussed again here 
in our new reversed perspective. The two muon beams in the close and far detectors and 
the oscillated electron beam at the far detector having widely different energy distributions, 
the effect of the reconstruction is expected to differ in all three. The muon neutrino energy 
distribution in the close detector, normalized with an energy integrated value of unity, 
$j, {E v ) is represented in Fig. [2] as a function of E v . At the arrival in the far detector it 
is reduced by a large factor which depends on the oscillation parameters and its expression 
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FIG. 2: (Color online) T2K flux energy distributions in the near detector (ND) and far detector 



(FD) for tho sets of oscillation parameters according to Ref. 



11[ and Ref. 



160- 



is : 

2 / &m 2 32 L 



4 cos 2 ^13 sin 2 8 2 3 (1 — cos 2 #i 3 sin 2 6*23) sin 



< D (E»J. (5) 



We use this expression as a definition of the far detector flux. We have kept in this expression 



the influence of the angle #13, which is now measured [15|: sin 2 #13 = 0.024 ± 0.004. Its 
effect is not totally negligible and it partly fills the dip of the energy distribution in the 
far detector. The oscillated distribution is shown as well in Fig. [2] for the values of the 
parameters of Ref. [16J and also for the best fit values of T2K [llj where the effect of #13 
is ignored. The products o{E v ^)§ Vii {E v ^) which represent the distributions of muon events 
before reconstruction in the close and far detector are shown in Fig. [3j We now apply 
our smearing procedure to these distributions and we obtain the smeared curves also shown 
in Fig. [3j The salient features are the broadening effects. In the close detector there is 
clear low energy enhancement, an effect of the multinucleon component. In the far detector, 
where the unsmeared distribution displays a pronounced dip, the smeared one acquires a 
low energy tail and the middle hole is largely filled, an effect also largely due to the np- 
nh cross section. All these smearing effects can be described as a tendency to escape the 
regions of high fluxes when one goes from true to reconstructed energies. We remind that 




FIG. 3: (Color online) T2K distributions per neutrons of muon events before (dashed lines) and 



after (continuous lines 



of the T2K best fit 



III. 



reconstruction in the near and far detector, evaluated with the parameters 



the opposite transformation from reconstructed to true energy tends instead to concentrate 
the events in the regions of high flux 4] . The ratio of the smeared distributions in the far 
and near detector displayed in Fig. H]is far from the unsmeared one, i.e., from the oscillation 
disappearance factor. 

Notice that the displacement in the far detector of the low energy peak towards smaller 
values by the smearing can be, to some extent, simulated by a decrease of the mass value in 
the unsmeared situation. As an illustration we compare in Fig. [5] the smeared curve with 
the T2K best fit, which has a mass value Am| 2 = 2.65 10 -3 eV 2 and the unsmeared curve 



with the best fit value Am 2 32 = 2.43 lCT 3 eV 2 of Ref. 



16| . The equivalence only holds in 



the first peak region but it is likely to that the inclusion of the reconstruction effects the 
analysis of the T2K disappearance data will result in a slight increase of the mass value. 
This investigation will be the object of a future work, following the approach of HQ. 
We will see that in the MiniBooNE case the equivalence between the introduction of the 
reconstruction effect and a lowering of the mass value in the unsmeared situation can be 
total. Turning now to the electron events distribution we display it in Fig. [6] without 
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FIG. 4: (Color online) disappearance oscillation probability with (continuous line) and without 
(dashed line) smearing. 
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FIG. 5: (Color online) The smeared distribution, D rec (E u ), with the T2K two-flavor best fit 11 | 



and the unsmeared curve, o{E U)1 )<fr Vli {E Vll ), with the best fit value of Ref. 



16j ]. The respective 



mass values are Am^ = 2.65 10 3 eV 2 and Am^ = 2.43 10 3 eV 2 . 
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FIG. 6: (Color online) T2K oscillation electron events energy distributions before (dashed lines) 
and after (continuous lines) smearing. The experimental histogram is also shown. 



and with the reconstruction correction. We also show the actual experimental histogram 



191 ] of significant events (i.e. with the background subtraction). The theoretical curve has 



been normalized to the same total number of these events. Again here the reconstruction 
correction tends to make events leak outside the high flux region, in agreement with the 
observed trend. 



B. MiniBooNE 

We turn to the MiniBooNE data which are now available for neutrino and antineutrinos 
. In both cases an electron excess has been observed and interpreted in terms of oscillations 
induced by the existence of sterile neutrinos. In the 3+1 hypothesis of only one sterile 
neutrino the electron neutrino flux is given by an expression similar to the usual two-flavor 
expression for active neutrinos. Starting from a muon neutrino energy distribution <§> v (E v ), 
the electron one generated by the oscillations is given by the following expression: 

/Am 2 T\ 

®u e {E Ve ) = ^(EJ sin 2 (2M sin 2 {-^j , (6) 
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FIG. 7: (Color online) Electron neutrino energy distributions as a function of E Ve from the oscil- 
lation expression of Eq. © for four set of values of the oscillation parameters. 



where Am^ is the difference of the square mass of the sterile neutrino and m\ (or m|) 
Short baseline oscillation experimental results imply large values of Am^, in the eV 2 range 
Constraints on this parameter have been provided by a series of data (see for examp 



Ref. [20] for a review or Ref. 21] ) and more recently by the ICARUS experiment 22]. 
Here our aim is to explore the oscillation phenomenon taking into account the problems of 
energy reconstruction, which has not been done previously. The data give the distribution 
of electron events as a function of the reconstructed neutrino energy. They have revealed a 
striking feature, denoted as the MiniBooNE anomaly [7, 8], namely an excess of events at 
low energies. There is instead a shortage of events above E Ue > 450 MeV. This has been 
the object of an intense debate. The low energy behavior of the data favors small values 
of the mass parameter, Am^ ~ 0.1 eV 2 , which concentrate the v e flux from the oscillation 
at low energies. But small values imply, in order to have enough events, large values of the 
sin 2 (2#4i ) which are not compatible with the constraints from other sets of data [20]. Figure 
[7] illustrates the electron neutrino beam energy distribution as a function of the energy 



obtained from the oscillation expression of Eq. fl6]) for some set values of sin 2 (2#4i) and 
Am^. Two values of the chosen set are reference values for LSND. We have used as the 
input of Eq. ([6]) the MiniBooNE muon flux <& v (E u ), as given in normalized with an 
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FIG. 8: (Color online) The electron distribution events o Ue {E Ve )Q Ve (E Ve ) corresponding to the 
fluxes of Fig. [71 



energy integrated value of unity. For large values of the Am^ parameter, the corresponding 
sin 2 term in the expression of the flux can take zero values in the relevant neutrino energy 
range, 0.2 GeV < E Ue < 1 GeV. For small values instead, it is a regularly decreasing function 
of the energy in the energy region above 0.2 GeV. When Am^ reaches values smaller than 
~ 0.1 eV 2 the sin function can be expanded and the energy distribution goes as 1/E% e . Once 
this situation is attained, for the concentration of the events at low energy there is no gain 
in lowering further the value of the mass, the energy distribution of the beam remaining the 
same, and it is possible to adjust the associated value of sin 2 (26*4!) to reach the required 
magnitude. We remind that the smearing effect does not depend on the value of sin 2 (26 , 41 ). 

The theoretical distributions of the electron events, the product of the "quasielastic" cross 
section by the oscillation electron neutrino energy distributions are displayed in Fig. [8]for the 
same oscillation parameters as Fig. [7J We now perform the smearing procedure according 
to the Eq. (j3J). Figure [9] illustrates its effect for a mass parameter Am^ = 0.45 eV 2 . The 
smeared curve is shifted at lower energies, the displacement of the peak is appreciable, ~ 100 
MeV. The smearing effect increases the strength at low energies and decreases it at high 
ones. This goes in the right direction, as too much strength is observed at low energies 
and not enough at large ones. It is important to stress that the smearing effect is here 
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FIG. 9: (Color online) Effect of the smearing procedure for a mass parameter Ara^ = 0.45 eV 2 
and comparison with the effect of the lowering, in the unsmeared case, of the mass to the value 
Am 2 = 0.1 eV 2 . 



equivalent and amounts to a lowering of the mass parameter. The smeared curve of Fig. 
can be reproduced in the unsmeared case with a value of the mass Am 2 x = 0.1 eV 2 as 
is illustrated in Fig. We have here adjusted the associated values of sin 2 (26> 41 ) to reach 
similar magnitudes. The equivalence is here complete. This means that, taking into account 
the smearing, a large mass value Am 2 ^ = 0.45 eV 2 allows the same quality of fit of the 
data than is obtained in the unsmeared case with a much smaller mass Am 2 x = 0.1 eV 2 . 
Obviously there is an important gain for the compatibility with the existing constraints since 
the angle 841 can be smaller with a larger mass value. The cases of smaller mass values, 
such as Am| x = 0.2 eV 2 , is also interesting. The corresponding smeared and unsmeared 
curves are displayed in Fig. [TUl showing the low energy shift, the peak position is lowered by 
~ 80 MeV. What is interesting is that it is impossible to reproduce, in the region of interest 
E v > 0.2 GeV, the smeared curve with an unsmeared one even taking a very small mass. 
The reason is that, as mentioned before, there is a limiting energy shape of the unsmeared 
curve, as is apparent in Fig. [JTjJ when the mass value is sufficiently small to reach the 
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and comparison with the lowering of the mass in the unsmeared case. 



1 / E% behavior for the energy dependence of the oscillated beam. The depopulation by the 
smearing of the energy region above 400 MeV at the benefit of the region below, even if it is 
not spectacular, and beyond what is achievable by a lowering of the value of the oscillation 
mass, is susceptible to improve the y 2 value of the best fit. The influence of the smearing 
effect on the data, which goes in the right direction to account for the low energy anomaly, 
is the main message of our study of the MiniBooNE data and it should have an influence in 
the analysis of the data. 

This being said, there is an additional problem which naturally comes to mind: the 
MiniBooNE muon neutrino beam contains a background of electron neutrinos which also 
undergo quasielastic events. These have to be evaluated and subtracted from the observed 
total number of electron events in order to extract the oscillation ones. These background 
electron events should be subject to the same reconstruction effects than the ones of the 
oscillated beam. This problem will be discussed in the following under the assumption that 
this background is evaluated in the same way as the oscillation events, namely from the 
knowledge of the electron neutrino background energy distribution and the "quasielastic" 
v e cross section. We are aware however that this treatment does not apply as such to the 
MiniBooNE results as the background evaluated in the actual analysis has been the subject 
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FIG. 11: (Color online) MiniBooNE electron events distribution for v e background in the unsmeared 
and smeared case. 



of a more detailed treatment than the one sketched here, with the information from the 
muon events. However our discussion may focus the attention on some of the problems 
which can occur in the analysis. 



1. Inclusion of the MiniBooNE v e background 

We want to investigate the correction introduced by the transformation from real to 
reconstructed energy in the background electron events distribution. In our previous work 
we have stressed the importance of the energy distribution in the reconstruction procedure. 
For the oscillation distribution, depending on the oscillation parameter, it is rather narrow. 
While instead the v e background has a broad distribution. As this distribution is different 
from the oscillation one, the reconstructed energy correction should be treated separately 
for the background and for the signal. We proceed for the background as for the signal. 
We start from the theoretical distribution for the u e background, cr(£' l/ J$ bacfc9rounrf (£' iye ). 



We have scaled here the background flux given in Ref. [23] by the same factor as the 
muon one, keeping in this way the right relative proportion between the background v e 
and the beams. We transform the (j(E Uf ,)^ background (E Ue ) in order to express it in terms 
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FIG. 12: (Color online) MiniBooNE muon events distribution in the unsmeared and smeared case. 

of reconstructed energy. The outcome is shown in Fig. [11] together with the theoretical 
distribution. As for the signal the smearing procedure increases the strength in the low 
energy region and reduces it beyond the peak position. This is essentially an effect of the 
np-nh piece, which is largely increased at small energies by the reconstruction procedure. We 
can therefore conclude that, in this description where the distribution of electron events from 
the v e background is evaluated theoretically as a product of the v e cross section and of the 
background flux, which means that the reconstructed energy corrections for the background 
are ignored , the electron events background is underestimated for low reconstructed neutrino 
energies E Ve < 0.6 GeV and overestimated for larger ones. Accordingly, in this picture, the 
oscillation signal excess obtained by subtraction of this background from the total signal is 
overestimated in the low energy region and underestimated in the high energy one, which 
would be of great interest. However this conclusion is entirely linked to the assumption made 
on the way in which the evaluation of the electron background is performed. Moreover, 
as this evaluation needs the absolute values, the cross section introduced is an essential 
ingredient. Here we have introduced the RPA cross section which introduces a quenching at 
low energies. It is partly compensated by the np-nh component, but the value of the cross 
section in the low energy domain remains an open question. 

In addition, as already mentioned, the actual evaluation of the electron background relies 



18 




FIG. 13: (Color online) Free (unsmeared case) and effective (smeared case) u- 12 C cross section per 
neutron for the MiniBooNE muonic events and for the background electron ones. 



in fact on a comparison with the muonic events, which also undergo the reconstruction 
correction. We give for completeness the reconstructed energy distribution also for the 
muonic events. The beam energy distribution is wider than the background v e one. The 
effect of the smearing is then even more pronounced, particularly in the low energy region 
where the distribution is appreciably enhanced, as is illustrated in the Fig. [T2j 

2. Effective cross sections 

All the smeared distributions D rec (E v ) can be expressed in terms of an effective "quasielas- 
tic" cross section afff(E v ) (which also includes the multinucleon emission) with the following 
definition, which applies to muons as well as to electrons: D rec (E u ) = a e J s (E V )$(E U ). We 
can now suppress the specification E v for the variable and write instead E v . From the eval- 
uation of the smeared distribution D rec (E v ) one deduces an effective cross section, which is 
different from the free one. This is the cross section that should be associated to the flux 
distribution in a theoretical evaluation to be directly compared to the data when these are 
expressed with the reconstructed energy. Its use incorporates the reconstruction correction. 
We stress once more that this effective cross section is not universal but it depends on the 
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FIG. 14: (Color online) Free (unsmeared case) and effective (smeared case) v- 12 C cross section per 
neutron for the MiniBooNE muonic events. The experimental MiniBooNE result [2] is also shown. 



particular beam energy distribution. Therefore on the calculational side there is no gain 
with respect to the previous evaluation, as this cross section differs for each energy distri- 
bution. Applied to the electron background, the effective cross section is different than the 
muonic one because on the one hand the lepton produced (e) is lighter and on the other 
hand the energy distribution of which it depends is also different. Without reconstruction, 
the electron cross section surpasses somewhat the muonic one, particularly in the threshold 
region. For the effective cross section the trend is reversed, the muonic one is larger than 
the background electron one, up to E v ~ 0.6 GeV. This difference of trends which emerges 
in this analysis shows some potential problems which may arise in the evaluation of the 
electron background from the contamination of the beam. In Fig. [THwe have also shown 
the MiniBooNE data points |2| for the muonic cross section which is expressed with the 
reconstructed energy. The points should then be compared to the effective cross section for 
muons. As compared to our first work |l) where the reconstruction correction was ignored 
the agreement remains satisfactory as expected, since it reflects the detailed agreement that 



we have js| for the double differential cross section . 
3. Antineutrinos 

Similar effects, although less pronounced, are present for antineutrino. They are not 
marked as in the neutrino case as, in our description, the np-nh events affect only the spin 
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FIG. 15: (Color online) Effect of the smearing procedure for antineutrino for a mass parameter 
Am]] = 0.45 eV 2 and comparison with the lowering of the mass in the unsmeared case to the value 
Am?i= 0.35 eV 2 . 



isospin response and are thus less influential in the reconstruction problem. For instance 
the smeared oscillation curve for a mass parameter Am| x = 0.45 eV 2 is equivalent in shape 
to an unsmeared one with Am 2 1 =0.35 eV 2 , as shown in Fig. [15j The gain is less than for 
neutrinos. But this awaits the test from the measurements of the antineutrino cross sections. 



IV. CONCLUSION 



We have discussed the application of our reconstruction method to several oscillation 
experiments. Starting from a theoretical energy distribution, product of the neutrino cross 
section by the neutrino energy distribution of the beam, we perform the reconstructed energy 
corrections. Their influence is such that the events tend to escape from the region of high 
fluxes with a tendency to concentrate at lower energies. The smeared distribution effect 
depends on the particular shape of the neutrino energy distribution, the correction being 
more pronounced for broad distributions. We apply our procedure to the three distributions 
measured in T2K: muonic distributions in the close and far detector and electron distribution 
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in the far detector. The effects are such that an analysis which takes into account the 
smearing effect is likely to lead to some increase of the oscillation mass value. 

We have also discussed the MiniBooNE results where the oscillations, if they exist, imply 
sterile neutrinos with a much larger mass parameter, in the eV range. The accumulation 
of electron events observed in the low energy sector favors relatively low values of this 
mass term which imply large mixing angles, not compatible with existing constraints. We 
have shown that this contradiction can be in part solved by the inclusion of the smearing 
effects. We have also pointed out some possible problems which may occur in the evaluation 
of the background from the v e contamination of the beam. In the antineutrino case 
where we predict similar effects but not as pronounced, the elucidation of the difference 
with neutrinos, which is important for the CP violation determination, will come from 
the measurement in progress of the double differential cross sections. In all instances the 
introduction of the smearing effect should improve the fit to the neutrino MiniBooNE data 
and its compatibility with the constraints from other sources. 
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